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Abstract 


Numerical  solutions  have  been  obtained  for  the  motion  of  a  gas 
bubble  in  an  incompressible  liquid  when  harmonic  pressure  oscillations 
are  imposed.  The  mean  pressure  in  the  liquid  was  taken  to  be  1  atm  and 
the  ratio  of  oscillating  pressure  amplitude.  A,  to  the  mean  pressure 

was  given  the  values  0.  2,  0.  5,  1.0  and  2.  0.  The  equilibrium  bubble 
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radius  was  chosen  to  be  Rq  =  10"  ,  10  and,  10"  cm,  and  the  angular 

frequency  of  the  pressure  variations  was  w  =  0—2  X  105,  0.7X  105,  and 

l  '■ 

1 .-2-X— TO5* per  sec.  The  phenomena  of  bubble  "explosion"  or  "collapse" 
were  found  for  large  A,  i.  e. ,  1. 0  and  2.  0.  For  the  smaller  values  of 
A,  the  bubble  radius  varied  with  time  within  well-defined  limits,  but 
nonlinear  effects  were  evident. 


Introduction 


The  equation  of  motion  for  a  gas  bubble  in  a  liquid  acted  upon  by 
a  variable  pressure  is  nonlinear,  and  solutions  for  the  radius  as  a 
function  of  time  cannot  be  given  explicitly  in  terms  of  simple,  or  familiar, 
functions.  The  example  of  a  pressure  variation  harmonic  in  time  is  ol 
particular  interest  since  it  corresponds  to  a  sound  field.  In  the  treat¬ 
ment  of  the  scattering  of  sound  by  a  gas  bubble  it  has  been  customary  to 
linearize  the  equation  of  motion  for  the  bubble  radius .  There  are,  how¬ 
ever,  many  experimental  situations  in  which  the  sound  pressure  ampli¬ 
tude  is  of  the  order  of  atmospheres.  The  acoustic  linearization  in  the 
liquid  remains  valid  with  high  precision  but  a  question  remains  regarding 
the  accuracy  of  the  linearization  of  the  bubble  equation. 

Numerical  integrations  of  the  bubble  radius  in  harmonic  pressure 
fields  have  been  carried  out  by  Noltingk  arid  Neppiras* .  The  primary 
difference  between  their  results  and  the  present  results  lies  in  the  great¬ 
er  range  of  parameters  used  here  and  also  in  the  longer  time  over  which 

the  bubble  histories  are  extended.  There  are  more  recent  calculations 

2 

by  Borotnikova  and  Soloukin  but  the  published  results  are  rather  limited 
and  more  complete  data  is  apparently  not  available  in  the  general  litera¬ 
ture.  Both  their  analysis  and  the  previous  work  of  Noltingk  and  Nep^i^as 
assume  that  the  bubble  oscillations  are  adiabatic.  It  is  known,  however, 

that  the  bubble  oscillations  in  the  range  of  concern  are  isothermal  rather 
3  4 

than  adiabatic  ’  and  the  calculations  presented  here  are  made  on  the 

5 

isothermal  basis.  Calculations  have  also  been  made  by  Flynn  but  his 
primary  concern  was  with  frequencies  much  higher  than  the  sonic  range 
which  will  be  considered  in  this  paper. 


I 

iw 


2 


There  are  two  physical  effects  which  contribute  elements  of 
indefiniteness  to  any  solutions  for  bubble  motions.  A  bubble  can  of 
course  be  set  into  free  oscillation;  in  the  limit  of  very  small  pressure 
amplitudes  such  oscillations  take  place  with  the  resonant  frequency  of 
the  bubble.  These  free  oscillations  are  excited  by  the  initial  conditions 
which  are  customarily  used  and  which  are  used  here.  The  bubble  is 
taken  to  be  in  equilibrium  and  at  rest  under  a  static  pressure  Pq  for 
t  <  0;  for  t>  0  the  oscillating  pressure  field  PqA  sinu>t  is  applied  and 
subsequent  motion  of  the  bubble  is  determined,  it  may  be  shown  explicitly 
from  the  linearized  solution  of  the  bubble  motion  that  viscosity  eventually 
damps  out  the  resonance  oscillations.  A  similar  behavior  is  to  be 
expected  for  the  nonlinear  solutions. 

No  consideration  will  be  given  here  to  possible  physical  mecha¬ 
nisms  which  give  a  stable  bubble  ot  "nuclei"  under  the  static  pressure 
PQ  for  the  indefinitely  long  period  t  <  0.  The  second  physical  effect 
which  is  not  considered  arises  from  rectified  diffusion  which  is  the 
process  by  means  of  which  a  pulsating  gas  bubble  grows  in  size  by 
transport  of  dissolved  gas  from  the  liquid  into  the  bubble.  This  rectified 
diffusion  is,  however,  a  slow  process  so  that  the  solutions  obtained  will 
remain  accurate  for  many  cycles  of  the  oscillating  pressure  field. 


Formulation  of  the  Problem 

When  the  motions  of  a  gas  bubble  are  sufficiently  slow  so  that  the 

effects  of  compressibility  may  be  neglected,  the  variations  of  the  bubble 

8  9 

radius  with  time  are  described  by  the  equation  * 


..  3  .  P^-Poo 

RR  +  4-  Rz=  - — 
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where  p  is  the  liquid  density,  p(R)  is  the  pressure  in  the  liquid  at  the 
bubble  wall,  and  p^  is  the  pressure  in  the  liquid  at  a  large  distance 
from  the  bubble.  It  is  assumed  that  effects  which  disturb  spherical  sym¬ 
metry  such  as  gravity  may  be  neglected.  The  condition  of  stress  contin¬ 
uity  across  the  bubble  boundary  gives 

p(R)  =  pg(R)  -  .  (2) 

In  Eq.  (2)  the  surface  tension  constant  is  c  and  p  (R)  is  the  gas  pres- 

S 

sure  in  the  bubble.  The  contribution  of  the  vapor  pressure  of  the  liquid 
to  the  pressure  within  the  bubble  has  been  omitted.  Since  the  bubble 
motions  will  be  isothermal  and  since  the  bubble  wall  velocities  are 
supposed  to  be  small,  the  vapor  pressure  may  be  taken  to  be  constant. 

A  constant  shift  in  the  pressure  level  is  not  significant  in  the  problem. 
The  viscous  stress,  -4{iR/R,  has  been  dropped  from  the  right  hand  side 
of  Eq.  (2).  This  term  has  a  significant  effect  on  the  bubble  motion  only 
after  a  long  time, as  has  been  remarked,  and  it  eventually  damps  out  the 
free  oscillations. 

For  t  <  0  the  bubble  is  at  rest  with  the  equilibrium  radius  Rq 

so  that 

p  {R  )  -  —  =  P  (3) 

®  o 


where  Pq  is  the  static  pressure  in  the  liquid.  The  equilibrium  relation 


calculations  which  have  been  carried  out  in  this  study  have  all  been  made 
with  Pq  =  1  atm.  For  t  >  0  the  pressure  at  a  distance  from  the  bubble 
is  supposed  to  be  disturbed  by  the  addition  of  a  harmonic  field  so  that 
the  total  applied  pressure  at  a  large  distance  from  the  bubble  is 
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pjt)  =  Pq(1+A  sinwt)  .  (4) 

The  bubble  motion  is  taken  to  have  the  initial  conditions  R(0)  =  R  , 

o 

R(0)  =  0.  The  gas  pressure  in  the  bubble  at  any  radius  is  determined  by 
the  isothermal  relation 

/  R 

Pg(R)  =  Pg(Rj[^-|  (5) 

Since  the  angular  frequencies,  o,  of  the  oscillating  pressure  will  lie 
in  the  sonic  range,  the  bubble  will  follow  the  isothermal  law  for  all  the 
radius  values  Rq  of  present  concern.  The  problem  is  completely 
defined  and  the  numerical  integration  can  proceed  in  a  straightforward 
way. 

It  is  convenient  to  introduce  the  dimensionless  radius 


(6) 

(7) 

(6) 

(9) 

(10) 


q(0)  =  0 
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Near  t  =  0,  -rj(t)  will  be  very  nearly  unity  so  that  for  small  time 
one  may  write 


i]{t)  =  1  +e(t) 


{12) 


where  |e(t)|  «  1.  Under  this  condition,  the  equation  obtained  from  (7)  for 
c(t)  may  be  linearized  and  integrated  analytically.  One  finds  readily  that, 
for  t  near  zero, 


T](t)  =  1 
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(15) 


is  the  isothermal  re  jnant  frequency.  This  initial  solution  (13)  is  used 
only  up  to  a  time  r  at  which  c  (t)  is  still  very  small.  The  solution 
which  is  developed  for  l  >  T  is,  of  course,  independent  ot  r  so  long  as 
T  is  sufficiently  small.  This  way  of  starting  the  numerical  solution  is 
for  convenience  only  and  serves  to  facilitate  the  development  of  the 
further  integration. 

Numerical  solutions  were  obtained  using  an  IBM  7094  Computer. 
The  numerical  integration  of  the  equations  was  performed  using  a  variable 
step  size  integration  scheme.  A  small  initial  step  size  was  given  as  an 
input;  the  integration  scheme  then  proceeded  to  integrate  each  step  with 
the  largest  step  size  possible,  consistent  with  the  specified  accuracy. 

It  was  found  that  different  choices  of  initial  step  size  had  no  effect  upon 
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the  character  of  the  solution. 

Discussion  of  the  Results 

The  results  are  presented  in  graphical  form  in  Figs.  1-9.  It 
may  be  observed  that,  as  might  be  expected,  the  largest  values  of  A, 
i.  e.  ,  1  and  2,  lead  to  large  variations  in  q  =  R/Rq  including  the  possibil¬ 
ity  of  collapse  or  explosion.  It  is  convenient  to  define  a  bubble  "collapse" 
as  a  change  in  tj  by  a  factor  of  10“  or  less;  similarly  a  bubble 
"explosion"  is  defined  as  a  change  in  r)  by  a  factor  of  102  or  more.  As 
a  general  result  it  appears  for  large  A  that  collapse  occurs  more 
frequently  than  explosion.  It  is  of  some  interest  to  observe  that,  for 
the  lower  frequencies  and  for  Rq  =  10  cm  and  10~*  cm,  no  collapse 
occurs  for  A  =  1.  The  rate  of  collapse  which  is  found  is  clearly  suf¬ 
ficiently  great  to  lead  to  cavitation  damage  and  cavitation  noise.  It  also 
is  evident  that  bubbles  with  radius  values  in  the  range  10  2  cm  to  10  4  cm 
can  act  as  nuclei  for  the  cavitation  phenomena  usually  observed.  There 
is  a  trend  for  the  higher  frequencies  to  give,  for  all  values  of  the  initial 
radius  and  of  A,  a  more  violent  radius  change;  there  is  a  trend  also  at 
higher  frequencies  for  collapse  or  explosion  to  occur  earlier  in  the 
bubble  history. 

The  bubble  radius  varies  within  well-defined  limits  for  A  =  0.  5, 
and  0.2.  In  the  figures  which  depict  these  solutions,  comparisons  with 
the  linearized  solutions  are  included;  where  they  are  omitted  the  two 
solutions  essentially  coincide.  For  the  high  frequencies  the  radius  varies 
in  a  more  erratic  manner  than  for  the  low  frequencies.  There  are  also 
indications  of  more  erratic  variations  in  radius  for  the  larger  initial 
radii.  The  resonance  oscillations  which  are  excited  have  the  effect  of 
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giving  a  nonperiodic  character  to  the  solutions. 

The  results  presented  are  for  a  phase  of  the  applied  oscillating 
pressure  which  gives  PqA  sinwt  with  A  >  0.  Essentially  the  same  behavior, 
as  is  shown  here,  is  obtained  for  A  <  0.  For  a  given  value  of  |  A j ,  the 
general  character  of  the  solutions  is  preserved.  This  behavior  is  to  be 
expected  since  the  effect  of  explosion,  or  collapse,  for  example,  appears 
only  after  a  few  cycles  of  the  applied  pressure.  The  specific  phase  at 
t  =  0  would  then  be  unimportant  provided  in  all  cases  R(0)  is  always 


small. 
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applied  at  t  =  0,  to  the  mean  pressure.  a>,  the  angular 

ncy  of  the  oscillating  pressure,  has  the  values 

10s  / sec  in  (a),  0,  7  X  105/sec  in  (b)  and  0.  2  X  105/sec 


Figure  5.  The  ratio,  of  the  bubble  radius  to  the  initial  radius,  R 

_3  C 

is  shown  as  a  function  of  time  for  R  =  10  cm.  The  mean 

o 

pressure  is  1  atm  and  A  is  0.5.  u>,  the  angular  frequency  of 

the  oscillating  pressure  is  1.2  X  105/sec  in  (a),  0.7X  105/sec 

in  (b),  and  0.  2  X  105/sec  in  (c).  The  dashed  curves  show  the 
linearized  solutions. 


I(U  10  «M) 


«t'07Xlcf  « 


Figure  6.  The  ratio,  rj(t),  of  the  bubble  radius  to  the  initial  radius,  Rq 
is  shown  as  a  function  of  time  for  R0  =  10“3  cm.  The  mean 
pressure  is  1  atm  and  A  is  0.2.  w,  the  angular  frequency  of 
the  oscillating  pressure  is  1.2X  itf/sec  in  (a),  0.7X  105/  sec 
in  (b),  and  0.  2X  lO^sec  in  (c). 
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Figure  8.  The  ratio,  t|(t),  of  the  bubble  radius  to  the  initial  radius,  R( 
is  shown  as  a  function  of  time  for  R0  =  10*4  cm.  The  mean 
pressure  is  1  atm  and  A  is  0.5.  w,  the  angular  frequency  of 

the  oscillating  pressure  is  1.2X  lC^/sec  in  (a),  0.  7  X  ltf/sec 
in  (b),  and  0.  2  X  1 0®  /  sec  in  (c). 
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Figure  9.  The  ratio,  t|(t),  of  the  bubble  radius  to  the  initial  radius,  Rq 
is  shown  as  a  function  of  time  for  R0  =  10^*  cm.  The  mean 
pressure  is  1  atm  and  A  is  0.2.  w,  the  angular  frequency  of 
the  oscillating  pressure  is  1.2X  lC^/sec  in  (a),  0.  7  X  10^/sec 
in  (b),  and  0.2X  itf/sec  in  (c). 
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